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ABSTRACT 

''  Much  current  practical  interest  exists  concerning  wave  propagation 
through  a  composite  medium  with  a  random  distribution  of  Inclusions:  parti¬ 
cles,  flakes,  long  continuous  or  chopped  fibers  in  a  homogeneous  matrix. 
Several  theoretical  studies  report  wave  speeds  and  attenuation  of  coherent 
plane  waves  propagating  through  an  elastic  homogeneous  medium  containing 
reinforcing  particles  or  fibers.  All  these  studies  assume  that  the  inter¬ 
face  between  the  matrix  and  the  inclusion  is  sharply  defined.  Also,  it  is 
mostly  assumed  chat  the  inclusion  is  perfectly  bonded  to  the  matrix. 

the  present  study  ,^*ej  analyzei/the  problem  of  damping  in  metal- 
matrix  composites  when  there  is  an  interface  layer  through  which  the  inclu¬ 
sion  property  changes  continuously  to  that  of  the  matrix. 

INTRODUCTION 

Determination  of  effective  elastic  moduli  and  damping  properties  of  a 
heterogeneous  material  by  using  elastic  waves  (propagating  or  standing)  is 
very  effective.  Several  theoretical  studies  show  that  for  long  wavelengths 
one  can  calculate  the  effective  wave  speeds  of  plane  longitudinal  and  shear 
waves  through  a  composite  material.  At  long  wavelengths  the  wave  speeds 
thus  calculated  are  non-disperslve  and  hence  provide  the  values  for  the 
static  effective  elastic  properties.  References  to  some  of  the  recent 
theoretical  and  experimental  works  can  be  found  in  [1-12].  The  scattering 
formulations  developed  in  [1-8]  provide  a  means  to  obtain  not  only  the 
effective  wave  speeds  but  also  the  damping  of  wave  amplitudes  due  to 
scattering. 

In  this  paper  we  present  results  of  some  of  our  recent  investigations 
of  phase  velocity  and  attenuation  of  plane  longitudinal  and  shear  waves  in 
a  medium  with  inclusions  with  interfacial  layers.  We  use  a  wave-scattering 
approach  together  with  uax's  quasicrystalline  approximation  to  predict  the 
phase  velocity  of  either  a  longitudinal  or  a  shear  wave  propagating  through 
the  medium.  The  scattering  approach  leads  also  to  an  estimation  of 
attenuation. 


1.  On  leave  from  Department  of  Mechanical  Engineering  II,  Tohoku  Univer¬ 
sity,  Sendai  980,  Japan. 
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MULTIPLE  SCATTERING  BY  A  DISTRIBUTION  OF  SPHERICAL  INCLUSIONS  AT  LOW  FRE- 
QUENCIES 


In  Chls  section  we  present  expressions  for  the  phase  velocities  and 
attenuation  of  longitudinal  and  shear  waves  in  a  medium  containing  a  dis¬ 
tribution  of  spherical  inclusions  with  thin  interface  layers  through  which 
the  elastic  properties  vary  rapidly,  but  continuously,  from  those  of  the 
inclusions  to  those  of  the  mattfx.  Such  interface  layers  are  often  present 
due  to  processing  [13,14]. 

In  order  to  use  the  scattering  approach  it  is  necessary  to  calculate 
the  scattered  field  due  to  a  single  spherical  inclusion  with  an  interface 
layer  (Fig.  1). 


Fig.  1.  Spherical  inclusion  with  interface  layer 

Consider  a  spherical  Inclusion  of  radius  a  and  elastic  properties,  A^,  p  , 
a  embedded  in  an  elastic  matrix  of  material  properties  Vi  -  and  P*. 


M.  eawaaeg  in  an  elastic  matrix  ui  material,  properties  ■'«*  ►*  and 
Also,  let  the  inclusion  be  separated  from  the  matrix  by  a  thin  layer  of 
uniform  thickness  (h(«a)  but  variable  properties  A(r),  M(r),  and  P(r). 
Here,  A,  v  denote  Lam4  constants  and  P  the  density.  Let  A(r),  u  (r)  be 
expressed  as 

A(r)  +  2p(r)  -  (A^  +  2^')  f(r),  a  <  r  <  a  +  h  (1) 


U (r)  -  g(r) ,  a  <  r  <  a  +  h 

f(r)  and  g(r)  are  general  functions  of  r.  A  special  case  arises  when 
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with  the  stipulation  that  f(r)  and  g(r)  with  their  first  derivatives  are 
continuous  in  (a,  a+h).  Since  h  is  assumed  to  be  much  smaller  than  a,  it 
follows  from  (3)  that  f'(a)  can  be  approximated  by 
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Note  that  X  '  and  U  '  are  Lam&  constants  of  the  interface  material  at  some 
value  of  r  (a<r<a+h;,  say  at  r-a+h/2. 

Another  special  case  would-be  that  the  Interface  material  possesses 
constant  properties.  Then  we  have 

f(r)  ■  g(r)  ■!,  a<r<a+h  •  (5) 
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We  also  make  the  assumption  that  h  is  very  much  smaller  than  the 
wavelength  of  the  propagating  wave.  Then,  to  first  order  in  h/X,  X  being 
the  wavelength. 
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Here  t . .  is  the  stress  tensor  and  superscripts  t,  s,  and  i  denote  the 
transmitted,  scattered,  and  incident  quantities,  respectively.  Note  that 
T8  ,  T..,  and  tT ,  appearing  above  are  calculated  at  r*a.  The  spherical  polar 
coordinates  r,  s,  <p  are  defined  in  Fig.  1.  Boundary  conditions  (6)  express 
the  fact  that,  to  first  order  in  h/X,  the  traction  components  do  not  suffer 
any  jump  across  the  layer.  However,  the  displacement  components  suffer 
jumps.  Two  parameters  that  characterize  interface  are: 

K  -  }  dx  K  -  /  -  —  (7) 

1  o  f<a+hx>  *  2  o  «(a+hx> 

Using  equations  (3)  and  (4)  in  (7),  it  can  be  shown  that  and  K2  are 
given  approximately  by 
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On  the  other  hand,  if  eq.  (5)  is  used,  then 

Kj  -  K2  -  1  (10) 

Mai  and  Bose  [2]  studied  a  problem  similar  to  the  one  considered  here. 
They  assumed  a  thin  viscous  third  layer  between  the  sphere  and  the  matrix 
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and  Imposed  the  condition  of  continuity  of  radial  displacement 


The  incident  wave  will  be  assumed  to  be  either  a  plane  longitudinal 
wave  propagating  in  the  positive  z-direction  or  a  plane  shear  wave  polar-* 
ized  in  the  x-direction  and  propagating  in  the  positive  z-direction.  Thus, 
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juVc.  and  k.*^/ c«.  u  denotes  the  circular  frequency  of  the  wave 
l’c2  l°a8itudlnal  and  shear  wave  speeds  in  the  matrix. 

The  factor  e  has  been  suppressed. 
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The  scattered  and  transmitted  fields  can  be  written  as 
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where  the  prime  denotes  that  k  and  k  are  to  be  replaced  by  k  '  («u)/c ')» 
and  k  '  («u)/c?' ) ,  respectively,  c,  and  c_'  are  the  wave  speeds  in1  the 
indu&on.  '  and  are  defined  in  [15]. 

The  constants  A,  B,  C,  A',  B' ,  C'  are  found  by  the  use  of  conditions 
(6).  For  details,  see  [15]. 

Once  A  ,  C  ,  and  B  are  determined,  the  scattered  field  is  then 
found  from  IquatTon  (12).  since  the  expressions  for  the  field  inside  the 
inclusion  will  not  be  needed  for  the  derivation  of  the  dispersion  equation 
governing  the  effective  wave  number  of  plane-wave  propagation  through  the 
composite  medium,  we  do  not  give  these  here.  In  the  following  we  present 
equations  governing  propagation  of  effective  plane  waves  through  a  medium 
composed  of  a  random  homogeneous  distribution  of  identical  spherical  inclu¬ 
sions  surrounded  by  the  layers  as  discussed  above. 

To  derive  approximately  the  phase  velocities  of  plane  waves  moving 
through  the  composite  medium,  we  assume  that  wavelengths  are  long  compared 
to  the  radius  of  each  inclusion.  In  this  long-wavelength  limit  it  can  be 
shown  that,  correct  to  0(eJ), 
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Expressions  for  P  ,  Q  ,  R  ,  and  S  are  defined  in  [15],  where  $  and  X 
•re  also  given.  11  a  a  Q  mn  mn 


Once  the  scattered  field  due  to  a  single  Inclusion  is  known,  multiple 
scattering  due  to  a  number  of  inclusions  can  easily  be  calculated.  In  par¬ 
ticular,  following  the  steps  discussed  before  by  us  [8]  it  can  be  shown 
that  effective  speeds  of  propagation  of  plane  longitudinal  and  shear  waves 
•re  given  by 


k*2  (1  +  9c?x)  (1  +  3cPq)  {1  +Y?2  (2  +  3t2)  } 

>  m  —  I.—  ■  ■  ■  ■  ■  —  . 

k2  1  -  15cP2  (1  +  3cPq)  +  |cP2  (2  +  3t2) 

k*2  (1  +  9cPj)  {1  +  |cP2  (2  +  3  t  2)  } 

k2  1  +  |cP2  (4  -  9t2) 

*  ,  *  *  * 

kl  "  ^Cl  *  k2  "  “/c2 

*  * 

c.  and  c  are  the  effective  wave  speeds  of  plane  longitudinal  and  shear 
waves,  respectively,  c  is  the  volume  concentration  of  inclusions  in  the 
matrix* 
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the  attenuation  caused  by  scattering  (to  this  order  of  approximation) 

can  also  be  calculated  using  equations  (12)  and  (15).  We  find  that  the 

attenuation  coefficients  ct  and -ct  are  given  by 
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It  may  be  noted  that  good  approximations  to  the  attenuation  coeffi¬ 
cients  are  obtained  if  P  ,  P.  and  P.  are  calculated  after  replacing  the 
matrix  properties  by  the  effective  composite  properties  given  by  Eqs.  (16) 
and  (17). 

SCATTERING  CROSS  SECTIONS  OF  A  SPHERICAL  INCLUSION  WITH  INTERFACE  LAYER  AT 
ARBITRARY  FREQUENCY 

Since  attenuation  coefficients  are  directly  related  to  the  scattering 
cross  sections,  in  this  section  we  present  some  results  for  scattering 
cross  sections  at  finite  frequencies  of  a  sphere  with  Interface  layer. 

From  equation  (12)  the  scattered  field  at  a  large  distance  from  the 
sphere  is  given  by 

u(8)  -  uP  +  u8  (20) 
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where  superscripts  p  and  s  denote  longitudinal  and  shear  wave  components, 
respectively.  It  is  easily  shown  that 
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where,  for  incident  longitudinal  wave  corresponding  to  the  first  term  in 
equation  (11), 
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The  scattering  cross  section  is  then 


8(0)  »  h~ 

00 

r 

"1 

n-3 

h(8)  - 

00 

z 

n-0 

t.  ?  <-i)n*L ) 


1 


lk 


1  n-0 


on 


For  incident  shear  wave  given  by  the  second  tens  in  Eq.  (11)  we  get 
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In  writing  (25)  we  have  rewritten  (16)  as 
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Here  subscripts  e  and  o  refer  to  even  and  odd  spherical  wave  functions  (see 
[16]). 


The  coefficients  appearing  in  (23)  and  (26)  were  calculated  when  the 
interface  layer  was  very  small  (h/a«l),  and  K.  and  K.  were  given  by  Eqs. 
(8)  and  (9).  Scattering  cross  sections  then  were  calculated  by  using  equa¬ 
tions  (26)  and  (27).  These  results  are  discussed  in  the  next  section. 


NUMERICAL  RESULTS  AND  DISCUSSION 

Attenuation  coefficients  for  two  composite  materials  were  calculated 
using  the  procedure  outlined  in  the  previous  two  sections.  The  first  com¬ 
posite  was  a  S1C/A1  material  which  was  studied  in  detail  in  [8].  The  con¬ 
stituent  properties  are  given  in  Table  1.  The  other  material  was  a 
lead/epoxy  composite  discussed  in  [17]*  Table  1  lists  the  properties  of 
the  constituents. 


Table  1.  Constituent  Properties  of  SiC/Al  and 
Lead/Epoxy  Composites 


Constituents 

P  3 
(g/cia  ) 

U 

(GPa) 

A+2jj 

(GPa) 

SiC 

3.181 

188.1 

474.2 

A1 

2.706 

26.7 

110.5 

Lead 

11.300 

8.35 

55.46 

Epoxy 

1.202 

1.71 

8.36 

Figures  2-5  show  the  variations  a  /k.  and  ctg/k  with  volume  concentra¬ 
tion  of  inclusions  at  different  frequencies.  Also^  shown  in  these  figures 
is  the  effect  of  a  thin  interface  layer  with  K.  and  K.  given  by  Eqs.  (8) 
end  (9).  For  SiC/Al  composite  Figs.  2  and  3  snow  that  the  effect  of  the 
interface  is  to  decrease  the  attenuation.  This  is  perhaps  understandable 
because  the  interface  is  not  a  sharp  discontinuity.  However,  as  shown  in 
the  case  of  Lead/Epoxy  composite  this  is  not  always  the  case.  These  fig¬ 
ures  also  show  that  attenuation  reaches  a  maximum  at  higher  and  higher  con¬ 
centrations  as  the  frequency  Increases.  Behavior  of  attenuation  vs. 


Attenuation  of  P  wave  in  SiC/Al  Composite  Fig.  3.  Attenuation  of  S  wave  in  SiC/Al  Composite 


Attenuation  of  P  wave  in  Lead/Epoxy  Composite  Fig.  5.  Attenuation  of  S  wave  in  Lead/Epoxy  Composite 
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concentration  in  Lead/Epoxy  composite  is  markedly  different  from  the  SiC/Al 
material.  It  is  seen  that  the  interface  generally  increases  the  attenua¬ 
tion.  The  peak  attenuation  is  reached  at  concentrations  that  first 
increase  with  frequency,  but  then  decrease  and  increase  again  at  high  fre¬ 
quencies.  Figure  6  shows  the  P-wave  attenuation  vs.  frequency  at  dif¬ 
ferent  concentrations.  This  agrees  well  with  the  results  obtained  in  [7] 
using  a  different  approach. 


,?S 


In  conclusion,  it  may  be  stated  that  attenuation  in  a  particulate  com¬ 
posite  depends  crucially  on  the  constituents  and  the  nature  of  interface 
layers.  These  layers  can  enhance  or  diminish  attenuation.  More  model  stu¬ 
dies  are  needed  to  characterize  attenuation  in  composites  with  thick  inter¬ 
face  layers.  This  work  is  in  progress  and  will  be  reported  later. 
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Fig.  6.  Attenuation  of  P  wave  in  Lead/Epoxy  Composite. 


